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I. INTRODUCTION 



Characterizing and quantifying multipartite entanglement is a problem whose complexity rapidly increases with the 
number of particles, and a major challenge within current quantum information science. In spite of intensive effort, a 
£C) ' complete understanding of entanglement properties remains limited to date to few-body small-dimensional composite 
\ quantum systems: in particular, such understanding has been achieved for pure states of three two-level systems 
(qubits) [H, 0, mixed-state entanglement having also been investigated for this system in Thus, the analysis of 
pure-state entanglement in an ensemble of four qubits is a critical test for any entanglement theory, as it is provides 
the first highly non-trivial case whose complexity remains tractable. Different approaches have been attempted so far 
for unraveling the classification of multipartite entanglement, ranging from so-called hyper-determinants 0], to normal 
forms Q, and invariants and covariants 0] of the relevant group of local transformations. Even if such methods 

offer equivalent answers for ensembles of two or three qubits, a complete description of four-qubit entanglement has 
p5 | only been obtained by Verstraete et al based on the method of normal forms, which simplifies considerably in this 
case thanks to the fact that the group SO(4, C) is isomorphic to SL(2, C) x SL(2, C). The resulting classification has 
been partially independently verified in [J]. Closely related with the problem of classification is, in turn, the problem 
of quantifying entanglement through appropriate measures, as the identification of proper classes should provide the 
physical boundaries for prospective measures. In addition, the invariants which are often utilized to discriminate 
among different entanglement classes satisfy themselves the minimum set of requirements that measures are expected 

CT 1 - to fulfill [H EH. 

In this work, we tackle the problem of pure-state four-qubit entanglement via a recently introduced approach 
based on nilpotent polynomials [Hj]. In addition to providing a simple entanglement criterion for any bipartition 
of a multipartite ensemble, the nilpotent method has the advantage of offering, in principle, a physically motivated 
procedure for entanglement classification, based on the idea of reducing the nilpotent polynomials to suitable canonic 
forms, which are extensive with respect to the number of subsystems and invariant under the desired groups of 
transformations. Such a reduction procedure is considerably facilitated if the dynamical equations of the polynomials 
are derived and employed. The coefficients of the resulting invariant forms have the same values as polynomial 
invariants, and may then be used for constructing measures of entanglement. 

The content of the paper is organized as follows. After recalling in Sec. [IT] the basic ingredients of the general 
nilpotent formalism, we specialize it in Sec. Illll to the four-qubit setting, and derive both gen eral and special entangle- 
ment classes for this ensemble. Note that we obtain more entanglement classes than in [lj|, as a consequence of the 
fact that we consider at each stage of our reduction procedure transformations that preserve the canonic form of the 
nilpotent polynomials. In Sec. IIV1 the problem of entanglement quantification is discussed in terms of the invariant 
coefficients of the nilpotent polynomials. Measures for comparing entanglement within classes are proposed, as well 
as a measure of genuine fourpartite entanglement. Sec. IV1 concludes with a summary of the results, and a discussion 
of the main advantages and limitations of our approach. 
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II. NILPOTENT POLYNOMIALS FOR ENTANGLEMENT DESCRIPTION 



Consider a pure state |^) describing an ensemble of n qubits. With respect to the computational basis in H ~ 
(C 2 )®", \iff) may be expressed in the form 

I*) = E {fei}=0 l^ fc " fc — 1 '" fcl \ k n k n-l ■ ■ ■ h) 

= Voo-o |00 • ■ -0) + ^io... |10-- -0) + . .. + V11-1 |H- ■■!), (1) 
where ipk n k n -i---ki S C. By introducing pseudospin creation operators cx^~, the above expression may be rewritten as 
|*) = Voo-o |00 • • • 0} + ip w ... a+ |0 • • • 0) + . . . + Tjj n ... 1 a+a+_ 1 ■ ■ ■ af |00 • • • 0) 

n s 

E ^-i-fall^) i°° • ■ -°> ' 

{fei}=0,l i=l ' 

that is, a polynomial in the nilpotent operator {erf} (recall that (erf) 2 = due to Pauli algebra), acting on the 
vacuum (or reference) state |0) = |00 • • - 0). By setting the population of the latter to be maximal (equal to one), we 
construct, equivalently, the nilpotent polynomial F, 

{k t }=0,l i=l {fci}=0,l W--- j=l 

Furthermore, by taking the logarithm of F, and by Taylor-expanding around the unit value of the vacuum state 
population, we obtain the nilpotential /, 

n 

f({af}) = ln[F({*f})]= • 

{fcj=0,l 4=1 

The nilpotential makes it possible to readily check whether two subsets A and B of qubits are entangled or not. The 
following criterion holds [12j |: 

The entanglement criterion: The subsets A and B of a binary partition of an assembly of n qubits are 
unentangled iff 

d 2 f{{ Xl }) 



dxkdx r 



, Vfc e A, Vm e B . 



Thus, A and B are disentangled iff /aub = + /b(Kb}). 

In spite of the fact that the nilpotential / gives the possibility of applying the above entanglement criterion, / 
may not yet be regarded as a satisfactory description of entanglement present in the overall composite system, as 
the latter should naturally be invariant under operations which act locally on individual subsystems only (see [Tij 
for a generalization of entanglement beyond the distinguishable subsystem framework we focus on here). The local 
transformations on each qubit may either be considered to be restricted to unitary transformations in SU(2) - in which 
case, we talk about su-entanglement - or they may be more generally allowed to be any invertible transformation 
in SL(2, C) - in which case, we talk about sZ-entanglement. Physically, the latter correspond to the family of 
stochastic local operations assisted by classical communication operations (SLOCC) [f| [l3j]. Under the action of 
local transformations (unitary or merely invertible), the state vector undergoes changes but still remains within a 
subset O, which coincides with a su-orbit (or, respectively, s/-orbit) within the overall Hilbert space TL. Thus, the 
nilpotential / should retain the same form for all states belonging to a given orbit, and a canonic form of the resulting 
nilpotential may accordingly be taken as an "orbit marker" . Canonic forms may be used as an alternative to the 
method of invariants [lH ] for identifying different orbits, thereby entanglement classes. The number of independent 
(real) parameters in a given canonic form should equal the number of independent invariants identifying the orbit, or 
else equal the dimension of the coset Ti/O. 

According to the general arguments given in [l2l . [l6j . the su-canonic nilpotential is defined as the nilpotential of 
the state in the orbit with the maximum reference state population. Under this condition, the orbit-marker is the 
canonic nilpotential, which we also term the su-tanglemeter, f c , 

fc(Wt}) - /W< + • • • , (2) 
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where the n linear terms are absent and the number of parameters involved equals the dimension of the coset, 
D su = 2 n+1 - 3n - 2, n > 3. 

In order to construct the sl-canonic nilpotential, or sZ-tanglemeter, we begin with the tanglemeter f c , and we further 
reduce the number of parameters down to D s i = 2 n+1 — 6n — 2, n > 4. To achieve this we impose the following 
conditions: in addition to the requirement for f c that all n terms linear in of be equal to zero, we require that all n 
terms of (n — l)-th order vanish as well. Thus, the sZ-tanglemeter, fc, takes the form 



Sifci^fl.n-l} 



(3) 



Since D s i < D su , different su-orbits may become equivalent under local Si-transformations. For this reason, the 

classification given by SL is more general than the one given by SU, thus usually the term "entanglement classes" is 
taken to refer to different sZ-orbits. 

Given an arbitrary pure state the task of determining the tanglemeter by applying local operations is, in 
general, not trivial. The difficulty is substantially reduced if one is able to take advantage of explicit dynamical 
equation obeyed by the nilpotential of the state, subject to appropriate consistency (or "feedback") conditions. For 
qubit systems, the dynamic equation reads 



i^Z = e -f He f 

dt 



where the generators of the local operations 



(4) 



(5) 



should be formally substituted as 



df 



*if = 



Of 



For the special case of local unitary operations, [P i ■) 
fc are 



dot 
(Pf) in Eq. 



where 



df 



dc 



(6) 



and the feedback conditions for obtaining 

(7) 

(8) 



are the coefficients of the linear terms in the nilpotential at a given time. 

A similar procedure for reducing the nilpotential to the canonic form fc may be carried out also for SL- 
transformations. We begin in this case by reducing / to the tanglemeter f c , so that the terms linear in &f vanish. 
Next, we apply SL operations as in Eq. ([5]), where however P~ and are no longer constrained to be complex 
conjugates, and choose such operations in such a way that the terms in the nilpotential involving the monomials of 
order one and of order n — 1 in af decrease exponentially with time. The two feedback conditions to be imposed in 
this case are: (i) the condition 



d 2 f 



dajdal 



(9) 



expressing P i via P^~ , which ensures that the nilpotential remains in the form of a tanglemeter at each stage; and 
(ii) the condition 



3 Wider" 



E P > 

m—1 



gn- 



df_ 

dot 



(10) 
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which ensures the exponential decrease of all n coefficients in front of the second-highest order terms. 

Unfortunately, no immediate physical meaning seems to be attributable in general to the requirements of vanishing 
of the s^-tanglemeter coefficients of (n — l)-th order - in contrast to the case of SU transformations, where vanishing 
of the first-order terms reflects maximum ground state population. Mathematically, however, such a requirement is 
suggested by symmetry considerations: n complex conditions are imposed on n complex coefficients of the same type. 
After having eliminated the monomials of orders 1 and (n — 1), it is possible to specify the scaling parameters Pf 
so that n additional conditions are imposed on the tanglemetcr coefficients. For example, we can set to unity the 
coefficients in front of the highest order term, and adjust (n — 1) coefficients in front of certain monomials to be equal 
to (n — 1) coefficients of other monomials. 

The condition in Eq. p0[) for P+ is written implicitly as a set of n linear equations that can be solved for generic 

states. However, no solution exists for those P + parameters corresponding to a zero determinant. Such singularities 
may correspond to special classes of entangled states which require separate consideration - as we are going to see 
explicitly in the four-qubit example. 



III. si-TANGLEMETERS FOR FOUR QUBITS 

A generic normalized pure state of four qubits may be described by 2 • 2 4 — 2 = 30 real parameters. The su- 
cntanglement of this state requires less parameters to be characterized, D su = 30 — 3 • 4 = 18 and, according to the 
discussion in Sec. UH for four qubits the su-tanglemeter defined in Eq. © reads 

f c = 3 <T+<7+ + f3 5 (7+(T+ + 9 CT+CJ+ +(3 6 (T+(J+ + 1O (7+(J+ +f3 12 (T+(J+ (11) 

+ 7 a^ at G i + Pi^t^t^i + 0ii (J X°' 2 (J i + 0u°~t a t a 2 + 0i^ a X G 3 a 2 a \ ■ 

In the above expression, we have used the local phase operations that did not contribute to the elimination of the 
linear coefficients to make the trilinear coefficients 7 , 0\ 3l fin, fin real numbers. In addition, a compact notation has 
been introduced by considering the indexes of fi as a binary representation of decimal numbers, e.g., 0011 1— » 3, etc. 

Allowing for more general local transformations on each qubit, such as indirect measurements with stochastic 
outcomes, the number of the parameters necessary to describe a state may be further reduced. The s^-tanglemeter 
|U of a generic state of four qubits contains D s i — 30 — 6 ■ 4 = 6 real parameters, and may be cast in the following 
form: 

fc = 03 Ka+ + a+a+) + fi 5 (a+o+ + a+a+) + fi 6 (a+af + afa+) 

+ (1 - fij - fi} - (3l)a+a+a+ai, (12) 

where the scaling factors (that is, the parameter in front of erf in Eq. (|3J)), have been chosen so that the fc becomes 
equivalent to the expression G a bcd in Theorem 2 of 0]. 

We proceed to explicitly illustrate the procedure for evaluating the sZ-tanglemeter in Eq. (fT2|) by means of the 
dynamic equations (J4j)- ((5j , starting from the su-tanglemeter given in Eq. (Ilip . First, one may notice that in the 
system of eleven first-order nonlinear differential equations for the coefficients (3i, the coupling of the second-order 
terms /3yC l +tT / *° * ne fourth-order term Px^a'l °~i occurs via the third-order terms fiia^oto-X ■> Pi^t a t a t > 
P11 o~t a 2 cr 2 ' Pu°~2 °~t a A ■ Thus, the time evolution of all /3$ stops when these third-order coefficients /?7, P13, (3, and 
/?14 vanish - indicating that for four qubits the sZ-tanglemeter is a stationary solution for the dynamic equations. If 
the coefficients P~ satisfy the requirement of Eq. (|9|) , which ensures that the nilpotential always remains in the form 
of a valid su-tanglemeter f c during such evolution, what it is left is to adjust the time dependence of the parameters 
P 2 + , P3 and P± so that they drive all four third-order coefficients to zero. 

From the differential equations of the third-order coefficients, 

ipu = -PfPxn + 2P+p 6 p w + 2P+p 6 p 12 + 2P+/? 10 /3 12 , 

i&a = 2P+l3 5 /3 9 -P+i3 15 + 2P+l3 5 p 12 + 2P+(3 9 p 12 , 

iPu = 2P+f3 3 f3 9 + 2P+p 3 f3 w -P+[3 15 + 2P+p 9 f3 w , 

ip 7 = 2P+f3 3 f3 5 + 2P+f3 3 f3 6 + 2P+p 5 p 6 -P+p 15 , (13) 

we see that, in the general case, feedback conditions may be imposed by a proper choice of the parameters P f , in 
such a way that these equations take the form 



$r = -Pr, ki = -0xi\ Pis = -l3i 3 ] $u = -0m 



(14) 
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The evolution implied by these equations brings, in turn, the nilpotential to the following form: 



(15) 



/ = f3 3 a+a+ + (3 5 (7+a+ + f3 9 a+a+ + (3 6 cr+a+ 
+ 0io44 + f3 12 cr^a^ + /3 15 a^a^a^a^ . 

We can invoke the four scaling operators e Bi(T i , and further reduce Eq. ([151) to the sZ-canonic form fc of Eq. (Ql 
unless one or more of the above (3 coefficients vanish. Such cases correspond to zero-measure manifolds - in other 
words, to special classes of entanglement. For example, when /3s = in (TT5|) . the tanglemeter may be cast, by scaling, 
in the form 



f (2) 

Jc 



44 + /%(44 + 44) 



+0 6 (44 + a+a+) + (1 - f3 2 5 - p 2 6 )a+a+<r+a+ , 
characterized by only two parameters. If f3 3 = 0io = 0, the sZ-tanglemeter reads 

fc ] = 44 + 44 + As(44 + 44) + (i - 6 2 )4444 

which only involves a single parameter. Lastly, if (3% — (3\q = f3g = 0, 



AO) 
Jc 



<7i (To Or, <7a 



(16) 



(17) 



(18) 



Note that the tanglemeters of Eqs. (fTrJ)) . (fTT)) . and (HHJ) correspond to the special families L a b C2 , L a2 b 2 and L a2 o 3@1 
of the classification given in Theorem 2 of [5J| . However, it is important to bear in mind that the latter classification 
applies to im-normalized states, whereas our tanglemeter corresponds to states of unit population in the reference 
state. 

When the fourth-order coefficient /?i5 = and, additionally, one or more of the quadratic coefficients are also 
zero, singular classes of states without genuine fourpartite entanglement emerge: for instance, the sZ-tanglemeter of a 
four-qubit W state, 

fc = 44 + 44 + 44> 

belongs to one of such classes, and separable states with tanglemeters of the type 

f c = 44 + 44 + 44 ! 

and similar, belong to other. 

On the other hand, reducing f c to the canonic sZ-form fc cannot be achieved when the determinant 



£> 4 = 



-/3i 5 2/? 6 /3io 2/3 6 /3i 2 2/?i /?i2 

20509 -/3i5 20 5 0i 2 20 9 12 

20 3 9 2/? 3 /3io -0i 5 20 9 0i O 

20 3 5 20 3 6 20 5 6 -015 



(19) 



of the system of differential equations (jT3J) vanishes - which makes it impossible to impose any required feedback 
conditions. In such a situation, we loose the functional independence of the right hand sides of (|13p . which ensures 
complete controllability of the dynamics of (3? , 0i 3 , 0n, and 0i4 in the generic case. In turn, this means that some 
linear combinations of these coefficients, determined by the system's eigenvectors, cannot be set to zero by any choice 
of , and a tanglemeter fc of a special form should be defined in such instances. In [l^], four special families of 
tanglemeters are derived, 



Jsi) 
Jc 



= (44 + 44) + 0s (44 + 44) 

+06 (afaf + 0-+0-+) 



T"3«2 U 1 - U 4 U 2 U 1 ~ ru A Ul i u l ~ 

+2 (0 5 6 - 3 6 + 3 5 ) 444^ 



444 



(20) 



As2) 

Jc 



+/? 6 (44 + 44) 



■07(444 -44°2 



-/3ii (44 4 



2444 <T i 



(21) 
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fc = a 2 a l + a 3 a t + a X a 2 + a i a 3 + a t a t + a 3 a 2 

+013 {^t (j t (j i + vt a i a t ~ a i a t a t ~ a i a t a t) 

+0n (cr 3 cr^cr+ - afa^t - v^a^a^ + cr 4~°3~ cr 2~) 

+2cr|cr+cr+CT 1 f . (22) 

/c = CT 3 _0 '2" cr i' + °4~ a t a t + a X a 2 a \ + aiat<*i 

+/3 3 a+a+ + /3 5 a+a+ + P 6 a+a+, (23) 

corresponding, respectively, to one, two, three, or four of the eigenvalues jj of T>4 vanishing, where explicitly 

71 = 0w - 2y/(3 5 f3 6 p g p w + 2^(3 3 (3 6 f3 9 Pi 2 - 2 v //3 3 /3 5 /3 1 o/3i2 , 

72 = 015 + 2^/(3 5 f3 6 p g p w - 2 v //3 3 /3 6 /3 9 /3 12 - 2 v //3 3 /3 5 /3 10 /3i2 , 

73 = As - 2V/3 5 /3 6 /3 9 /3io - 2V/3 3 /3 6 /?9/3i2 + 2V/3 3 /3 5 /3io/3i2 , 

74 = Pis + 2^/f3 5 f3 6 p g [3 w + 2 v //3 3 /3 6 /3 9/ 3 12 + 2^f3 3 f3 5 p w (3 12 . (24) 

Observe that the number of parameters in such special tanglemeters is 3 complex numbers, the same as in the general 
case given in Eq. (fT2"|) . 

At present, it remains to be proved whether all four special tanglemeters (f20 |) -(|23 |) correspond to distinct special en- 
tanglement classes, since they result from considering a dynamic evolution based on a series of sequential infinitesimal 
local operations which preserve the su-canonic form of the nilpotential. Thus, a situation where some of the obtained 
tanglemeters turn out to be equivalent under a finite local s^-transformation, cannot be ruled out in principle by our 
current approach. In [f|, such special classes are not explicitly identified, although the last three classes of Theorem 
2 - that is, io 7ffiI > ^(^iCW , and Lo 5lsS , may be easily identified as special cases of Eq. (|23|) when one or more terms 
vanish. The class L a \, 3 in [f| is not identified by our method. 

We summarize in Table 1 the entanglement classes for pure states of four qubits we have thus obtained. 



IV. ENTANGLEMENT MEASURES FOR FOUR QUBITS 
A. Measures for si- and su-entanglement 



From an information-theoretic standpoint, the construction of well-defined entanglement measures typically relies 
on the concept of entanglement monotone, that is, of a quantity that is required to be invariant under local unitary 
transformations and non-increasing on average under LOCC transformations (l3| . For instance, the most widely 
utilized measures for two and three qubits, the concurrence, C, and the residual entanglement (or 3-tangle), r, [171 ] 
are entanglement monotones. For a four-qubit system, we have seen in Sec. Mil that the classification is much richer 
than in the case of three qubits. In the context of such a classification, we would like to first revisit the role of 
entanglement monotones, and then argue that another class of measures may also be meaningful. In particular, 
we show how a measure for fourpartitc entanglement should be also more precisely defined by imposing additional 
requirements beside the ones mentioned above. 

A standard way to construct entanglement monotones is based on exploiting polynomial (algebraic) invariants. 
Polynomial invariants are polynomial functions of the state coefficients, and a linearly independent finite set of them 
may be used to distinguish different orbits in the same way the set of invariant tanglemeter's coefficients does. For 
example, for a three-qubit system, five (as many as the tanglemeter's parameters) independent invariants under local 
unitary transformations exist namely the three real numbers 

/ 3 = 4>r 3k r iJP ^ m n P r mnk , (25) 

and the real and the imaginary part of a complex number, 

h + i/ 5 - ^jk^^mnr^ ■ (26) 
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In the above equations, i\) l i k = e** e ipi'j'k', with the convention that summation over repeated indexes ranging 
over {0, 1} is left implicit, ip* 1 ^ denotes the complex conjugate of tpijk, and e M is the antisymmetric tensor of rank 
2. The sw-invariant quantity t = 2\I± + is exactly the 3-tangle, which also remains invariant under the class 
of local transformations (8)?_ 1 5£j(2, C). It was proved in [j| that sZ-invariants behave as entanglement monotones 
for normalized pure states, since the (square) vector length ^ ipiip* is non-increasing under sZ-transformations and 
thus may be employed as a measure of su-entanglement within a given si orbit. The main reason for choosing such 
a vector length as a measure is based on the relation between the determinant D et < 1 of the physical transforma- 
tion corresponding to the chosen sZ-transformation and the probability ipiip*)~ l of the desired outcome of the 
indirect measurement implementing this transformation: The (n — l)-th power of the probability upper-bounds the 
determinant, D et < CiV'iV'i) -71 - 

Note that, by definition, an entanglement monotone is an object able to quantify su-entanglement by distinguishing 
different sit-orbits that belong to the same sZ-orbit. However, in the case of four (or more) qubits, there exists an 
infinite number of general sZ-orbits (see Sec. IIII[) . This suggests that measures able to compare the si -entanglement 
between such general orbits should be considered in addition to the su-measures. A reasonable suggestion for sl- 
cntanglement measures is provided by sZ-invariants that are also scaling invariants and, therefore, are independent of 
the specific normalization of the state. One may construct sZ-invariants for a four-qubit ensemble in a way similar to 
how the invariant I4 + il$ of Eq. (|26[) is constructed; that is, by taking products of several factors ~ tf> (but not factors 
~ tp*) and by considering contractions over S'[/(2)-indexes with invariant antisymmetric tensors e M . The simplest 
combination one finds in this way, 

I {2) = V«W^ iJW , (27) 
is a sZ-invariant of second order. There also exist three different sZ-invariants of fourth order, 

Il4 =^3 = ^iklrt imn Homnp^° klP ■ (28) 

The ratios 1$ /(I™) 2 , 1$ /(/ (2) ) 2 , and 1$ /(I™) 2 are, in addition, invariant with respect to multiplication of the 
state vector by an arbitrary complex constant. Were these ratio linearly independent, they would suffice for a complete 
characterization of four-qubit entanglement. However, they are not. The following identity, 

= (29) 



makes such quantities inconvenient for entanglement characterization. 

Thus, it is necessary to turn to the sixth-order invariants. We consider the following three independent combinations, 

-^12 g {.^Pingd^Prarko^P sjph ^Pingo 

^23 ~ Q i^ijpo^mngh^srkd ' i Pijpd1pmngo' l Psrkh) 1p ^ 1p' ^ 1p ^ ' , 

Il3 =\ (iPvkh^mnpdlpsrgo - ^ijgh^mnkd^srpo)^™^^^ , (30) 

whose differences give the invariants of Eq. (|28p multiplied by l' 2 '. The explicit form of these invariants for a generic 
state is awkward. However, they take a simple form for the canonic state under sZ-transformations, which allows us 
to explicitly relate them to the canonic amplitudes. One finds 



^0000 - y2(/(2)) 1/4 : 



01100 = 00011 = 
"01001 = 00110 = 
■00101 = "01010 = 











2(/( 2 ) 






■y/iiV+Q- 


-^i[V+Q+(i^f /2 




2(/( 2 ) 


sl/4 






-V / ^ ) +Q+(/< 2 ») 3/2 



2(7(2))' 
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where Q is a root of the following cubic equation: 

(4f + Q)(4f + Q)(4f + Q) = (i {2) ) 3 Q 2 • (32) 



The set of Eqs. (|31[) determines the canonic state vector form with respect to pure SX-transformations. Upon 
dividing Eqs. ([31]) by V'oooo, the ratios V'noo/V'oooo, V'iooi/V'oooo, and V'oioi/V'oooo respectively yield the sZ-tanglemeter 
coefficients fa, fa, and fa, which are also scaling- invariant. Different roots of the cubic equation (|32l) yield different 
sZ-canonic states related by SL transformations. We can choose one particular root by minimizing the difference 
between the normalization of the canonic state and the initial normalization. Thus, as conjectured in Sec. IIII1 the 
sZ-entanglement in the four-qubit assembly may be completely characterized by three independent scale-invariant 
complex ratios, 



fa 



fl V / V / ^-v / ^-V / ^+(^) 3/2 

P6 = / ; ; ; , (33) 

emerging from the invariants of Eqs. (|30p - (|27[) . In view of this, a natural measure of sZ-entanglement is provided 
by the sum of squared moduli of the sZ-tanglemeter coefficients fa S2 — |/9| 2 - This yields S2 = for the GHZ 
canonic state, whereas S2 ^ for all other states, thereby exhibiting a similar behavior to the hyper-determinant (4j. 
Accordingly, this measure quantifies how close the orbit is to the GHZ-orbit. The quantity ^ \(3 — fa\ 2 may likewise 
serve as a measure characterizing the distance between two different sZ-orbits. 

As a next question, we wish to suggest a simple measure for characterizing su-entanglement in four qubits. A 
natural candidate is the sum S\ = ^ \ip | 2 over the probabilities in Eq. (|31[) . which gives the standard normalization 
of the canonic- like state. Once the invariants of Eqs. (|27p - (l3"01) are calculated for a state with unit normalization, this 
sum quantifies the extent by which the SL transformation required for setting the state to the canonic form differs 

from a unitary transformation. Thus, ln^ \ip\ 2 provides us with a suitable measure of such a non-unitarity. By 

construction, the latter quantity is able to discriminate between different su-orbits that belong to the same sZ-orbit. 
Interestingly, as found in (T^ |. this measure exhibits a strong correlation with the quadratic sZ-invariant I^ 2 \ 



B. Measures for fourpartite entanglement 

Having suggested si- and sit-measures for four qubits in terms of the tanglemeter's coefficients, we finally proceed 
to address the more delicate issue of constructing a measure of genuine fourpartite entanglement [18j . In addition to 
behaving as an entanglement monotone, such a measure should satisfy the requirement of being zero in the sZ-orbits 
that do not bear genuine fourpartite entan glem ent. Within constructions based on sZ-invariants (different approaches 
have also been suggested, see e.g. [l8l. [l9lT20j). the combination of invariants able to satisfy the last requirement is 
not known to date. For example, 1^ in Eq. ([2T)l is a low-order entanglement monotone, but it cannot serve as a 
good fourpartite measure since it attains its maximum value 1 for both the four-qubit GHZ state and for a product 
of two Bell pairs, that is a four-qubit state which manifestly contains no genuine fourpartite correlations [Hj. The 
so-called 4-concurrence introduced in [2l|, that is just |^ 2 '|, exhibits a similar unfavorable behavior. On the other 
hand, the hyper-determinant A p[ is nonzero in the general family of orbits G a bcd, and zero in all others as well as 
in the GHZ orbit. According to our results (Table I), recall that the families of orbits L a & C2 , L a2 b 2) and L a2 o mi are 
derived as special cases of the general family, and also contain genuine fourpartite entanglement in general. 

Observing that the determinant P4 of the infinitesimal transformations given in Eq. (|19p is precisely equal to zero 
in the orbits G a bcd, L a b C2 , L a2 b 2 , and L a2 o 3ei , we express it in terms of the canonic state amplitudes, 

*15 2lp 6 1p w 2-06^12 2^10^12 
2^ 5 V>9 *i 5 2i/j 5 1pi2 2l/j 9 1pi2 
2^ 3 V>9 2^3^10 #15 2^9^10 

2ip 3 ip 5 2ip 3 ip 6 2ip 5 ip 6 *i 5 



(34) 
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where ^15 = — V'isV'O + ^6^9 + ^z^n + V'sV'io- Our proposal is to consider the quantity 

Ki = VMM 4 (35) 

as a measure of proper fourpartite entanglement. Here, the normalization A = JT IV'i/'/'ol 2 ensures that JC4 con- 
sistently ranges between 0, 1. Note that /C4 is constructed as a function of su-canonic amplitudes, thus it remains 
invariant under local unitary transformations, while it is not a sZ-invariant. 

In order to gain further insight into the properties of /C4, its behavior on a typical four-qubit state is numerically 
investigated. The distribution of /C4 for pure states sampled uniformly over the Haar measure is depicted in Fig. [I] A 
non-monotonic behavior is clearly seen, peaked around the value 0.1, qualitatively resembling the distribution of the 
3-tangle for three-qubit pure states [12]. Unlike the 4-concurrence, /C4 attains its maximum value 1 for the GHZ state, 
and gives zero for all the states which are separable in some way. For a four qubit cluster state [23[, for instance, it 
gives K.4 = 0.3265, which is well above the average value. 
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FIG. 1: The distribution of K.4 for a set of 10 4 four-qubit pure states sampled uniformly over the Haar measure. See for 
comparison Fig. 2 (inset) in [22|. 

While the above features make IC4 an attractive candidate for quantifying fourpartite entanglement, a main disad- 
vantage of /C4 is that it inherits the redundancy of our classification, vanishing whenever the general class of orbits 
cannot be reached by infinitesimal transformations - irrespective of whether it might be reached by finite transfor- 
mations. Furthermore, the calculation of /C4 for a given pure state requires in general that the latter is first reduced 
to its su-canonic form. On the other hand, extending the construction of this measure to n > 4 qubits is relatively 
straightforward in principle. For n = 4, the fact that the determinant K4 does not contain the second-highest order 
terms is a sign that this measure is approximate for arbitrary states. However, this effect may expected to become 
less pronounced (hence the accuracy of such approximation improves) with increasing n. Finally, it remains an open 
question to ascertain whether this measure is an entanglement monotone. 



V. DISCUSSION 



In summary, we have demonstrated how the approach based on nilpotent polynomials may be employed to iden- 
tify entanglement classes for the illustrative yet highly nontrivial situation of four qubits in a pure state. Even if 
the approach is somewhat redundant compared to mathematically more sophisticated methods, we believe it has 
the advantage of offering a clearer physical interpretation, and it may also be extended straightforwardly to larger 
multipartite ensembles and higher-dimensional subsystems. 

In the context of the resulting classification, we have suggested additional class of measures beside the existing 
ones, which remain invariant under either local unitary (su) or arbitrary local invertible (si) transformations. We 
employ the nilpotent invariant coefficients for the construction of such measures as an alternative route to invoking 
polynomial algebraic invariants. 

Finally, we have also suggested a measure of genuine fourpartite entanglement. Our prospective measure is both, by 
construction, a su-invariant and it vanishes on the special orbits where no genuine fourpartite entanglement exists. We 
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also note that while it would be very illustrative to apply the nilpotent method in the well-explored three qubit case, 
the procedure implemented here is not viable, since in this case the SL coset dimension D s i < and, consequently, 
the dynamical set of equations employed in the analysis becomes degenerate. 

It is our hope that the results presented here may serve as a stimulus to prompt further investigations and appli- 
cations of the nilpotent polynomial formalism as a tool exploring entanglement. 
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General class 



3 complex parameters 



G a 

Singular 3D classes 



+fe(a+a+ + a+a+) + (1 - fe 2 - fe 2 - a+ a+ a+ 

3 complex parameters 



G b 



G d 



Singular 2D classes 



f = (3 :i (a+a+ + a+a+) + fe (a+a+ + a+a+) + fe (a+a+ + a+a+) 
+af at at - at at at + at at at - at at at 
+2 (fefe - fefe + fefe) at at at at 

f = °t a t + a t a t + a 2 G t + a t a 4 + fe ( a t a t + a t a t) 

+fe (at a 2 a t - a 2 a t a 4 ) + fei (at at at - a t a t a t) 
+2atatatat 

f = a t a t + a t a t + a 2 a t + a t a 4 + a t°t + a t a t 

+/3i4 (at atat - at atat + at atat - atat at) 
+/3i3 (at atat + at atat - at atat - atat at) 
+/?n (at atat - at atat - at atat + atat at) 
+2a+a+a+a+ 

f = atat at + at at at + at atat + at at at + 
feat at + feat at + feat at 

2 complex parameters 



LG2 a 
LG2 b 
LG2 C 
Singular ID classes 



+ folptet + + fc(°t<rt + °t°t) + (1 - 0s - Pl)atatatat 

f = atat + atat + fe{at a+ + atat) + fe{at a+ + atat) 

f = a t a t a t + a t a t a t + a t a t a 4 + a t a t + fe<*t a t + fe a t a t 
1 complex parameters 



LGl a 
LGl b 



f = atat + atat + fe(at at + atat) + (1 - /^K+a+a+o^ 

/ = a t a t a t + a t a t a 4 + G t°t + a t a t + fe°t°t 



point classes 






no parameters 


S a 




f-- 


= at °t a t + a t a t a t + °t a t 


s b 




/ = 


a t a 2 a t + a t a t a 4 + ot a t a t 


S c 






f = °~t a t a t + °t a t a 4 


s d 






c + + + 

/ = a 1 a 2 a A 


S e 


f 




at + atat + atat + at atat- 


Sf 


f 


=*t 


at + atat + atat + at at at- 



TABLE I: Classification of four-qubit entanglement classes following from SL(2, C) transformation properties of the canonic 
form, see Sec. 3. 



